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Abstract 

We study general zero range processes with different types of particles on a d- 
dimensional lattice with periodic boundary conditions. A necessary and sufficient 
condition on the jump rates for the existence of stationary product measures is 
established. For translation invariant jump rates we prove the hydrodynamic limit 
on the Euler scale using Yau's relative entropy method. The limit equation is 
a system of conservation laws, which are hyperbolic and have a globally convex 
entropy. We analyze this system in terms of entropy variables. In addition we 
obtain stationary density profiles for open boundaries. 

1 Introduction 

The zero range process is a stochastic particle system on the ci-dimensional lattice 7L d 
where the jump rate g(k) of a given particle depends only on the occupation number k 
at its current position. This model was originally introduced as a simple example of an 
interacting Markov process jl]. Various properties have been established, among them 
the existence of the dynamics under very general conditions, classification of invariant 
measures, and hydrodynamic limits [3 13 El • 

In this paper we generalize the zero range process to n different types of particles. The 
jump rate gt of the i-th component depends on the occupation numbers of all n species 
at a given site. For the process to have stationary product measures, these rates cannot 
be chosen arbitrarily. We find that such measures exist if and only if the logarithm of 
the jump rates is given as the lattice derivative VjG of a suitable potential G. In 
general, the stationary product measures do not factorize with respect to the different 
components. 

Our goal is the derivation of the hydrodynamical equations on the Euler scale, which 
are expected to be given by a system of n conservation laws. For this purpose we use 
the relative entropy method of Yau [3], which, while directly applicable to the present 
case, has the disadvantage to yield the desired result only up to the first shock. So far we 
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did not attempt to extend our result to all times following e.g. the lines in j^j. We show 
that the system of conservation laws is hyperbolic and the thermodynamic entropy of the 
stationary measure is a Lax entropy [7j. This property follows from certain reciprocity 
relations for the steady currents of the components. Such relations have recently been 
established for one-dimensional interacting particle systems restricted to nearest neighbor 
interaction and finite occupation numbers jH]. 

One source of motivation for this work comes from the analysis of one-dimensional 
driven diffusive systems with open boundaries. For one-component systems the analysis of 
the hydrodynamic limit equation led to the theory of boundary-induced phase transitions 
which provides a general framework for a quantitative description of the steady-state 
selection in systems which are in contact with particle reservoirs at their boundary [9, 
ITU] . In systems with more than one conserved quantity interesting new phenomena have 
been found such as phase separation and spontaneous symmetry breaking flj |Q, for 
a recent review see [13J. Again it is natural to ask for principles governing steady-state 
selection and the resulting phase diagram in systems with many species of particles. 
The macroscopic behavior of such systems has been examined to some extent fl] only 
recently and there are very few rigorous results (HI ■ I n our note we analyze the system 
of conservation laws for open boundaries in terms of entropy variables ^H]- For stationary 
solutions the system decouples, and we are able to derive stationary density profiles, which 
we state explicitly for one-dimensional geometry. 

2 Stationary product measures 
2.1 The n-component zero range process 

Let us consider a zero range process with translation invariant jump rates on the d- 
dimensional torus = (Z/LZ) d . There are n different species of particles and let 
r/j(x) G N be the number of particles of component i G {1, . . . , n} on site x G A L , where 
N = {0,1,2,...}. The state space is given by VLl = (N n ) AL and we denote a particle 

configuration by rj = ( r ?( x )) xgAi = ((?7i(x), • • • > ^(x)) j . At a given site x G A^, 

the number //j(x) of z-type particles decreases by one after an exponential waiting time 
with rate jjm(x)J and the leaving particle jumps to site x + y with probability Pi(y). 
The jump probabilities pi : Z — > [0, 1] are normalized and assumed to be of finite range 
R G Z+ i.e. 



To exclude hidden conservation laws the p^s have to be irreducible, so that every particle 
can reach any site of the lattice with positive probability. The rate function : N n — > 
[0, oo) vanishes for all k = (ki, . . . , k n ) G N n with ki = and is otherwise positive and 
uniformly bounded from below, i.e. using the shorthand k = (k 4 , ki), 




and 




(2.1) 



y 




and 
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for all i = 1, . . . ,n. With these assumptions the generator of the zero range process is 
given by 

n 

oc/)fa)= E E^(^( x ))^(y)(w ;x ' x+y )-/(^))' ( 2 - 3 ) 

x,yGA L i=l 

regarded as a linear operator on C($7l, R). The configuration r7* ;x,x+y results from 77 
after one particle of component i is moved from x to x + y, i.e. rjj* y (z) = 77^ (z) + 
^i(^z,(x+y)modL d ~~ ^z,x) for all z £ A L , j = 1, . . . ,?7,. The number of particles N^rj) = 
SxeAi^( x ) °f component i is conserved, and these are the only conserved quantities. 
They divide the configuration space into finite, invariant subsets with fixed iVj G Z + , 

2 = 1, . . . , n. Restricted to such a subset, £ is a finite dimensional matrix and the process 
is well defined. However, for L = oo this is true only for "reasonable" initial conditions 
and under additional assumptions on the jump rates, which are given in Section 3.1 (cf. 

2.2 Existence of stationary product measures 

For the one-component process n — 1, as well known, there exists a family of translation 
invariant stationary product measures (see e.g. [21 131 IT7j) 

1 HW .J oo l 

n W) ""» n ^ ^^)=E«"n^- ^ 

xeA L fc=l yiv ' 1=0 fc=l yiv 7 

The parameter \i e R is the chemical potential and controls the average particle density, 
the normalizing constant Z{ji) is the partition function. Often /i is replaced by the fu- 
gacity = exp[/i] 6 [0, oo). We use this notation in Section 4.2, but for Sections 2 and 

3 the chemical potential turns out to be more convenient. In the case n > 1, the sta- 
tionary measures are of product form only under the following condition on the jump rates: 

Assumption: For every i,j G {l,...,n} and k = (ki,...,k n ) G N n with ki,kj > 
let 

9l (k) 9j (k\ ki-l)= 9j (k) gi Qt>, kj - 1) . (2.5) 



This assumption is equivalent to the existence of a potential G : N n — > R for the 
logarithm of the jump rates such that 

logg i (k) = G{k)-G(k\k i -l) . (2.6) 

Given G, the jump rates defined via (|2.6j) clearly satisfy (|2.5|) by construction. On the 
other hand for given jump rates obeying (|2.5p one can define G recursively via (|2.6|) 
by fixing G(0, . . . , 0) =0. This construction does not depend on the order of summation, 
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since by (J2.5|) the sum over every closed path vanishes. For example one can choose 
ki 

G(k) = ^log£ 1 ( Jl ,0,...,0) + 

h=l 

k 2 k n 

E I°g92(h,j2,0, ...,0) + ... + ^ logg n (kx, ■ ■ -,k n -iJn) ■ (2-7) 

J2 = l jn = l 

Theorem 1: (Stationary product measures) 

The zero range process defined in \2.§) with more than one component, n > 1, has sta- 
tionary product measures if and only if the condition \2.5]) . equivalently \2.b)) . is fulfilled. 
In this case the family of stationary measures can be written as 



xeAi, 



8=1 



(2.8) 



with the chemical potentials /x = (/ii, . . . , fi n ) G as parameters. is the domain of 
convergence of the partition function 



Z{n) = ^ exp - G(k) + ^fi i k i 



(2.9) 



keN™ i=l 

Dn is a nonempty, convex subset ofW 1 with infinite volume measure. 



Proof: First we assume that ()2.5|) and ()2.6|) are satisfied and show that vf^ defined 
in ()2.8|) is stationary. Since G exists by assumption, the measure uf^ is well defined on 
and {/x G R n |/ij < hgg*, i = 1, . . . ,n} C because of For Z(fi l ), Z(n 2 ) < °o 

it is easy to see that Z(qn l + (1 — g)/x 2 ) < oo for all q G [0, 1], thus is convex. The 
remainder of the first part of the proof is a straightforward generalization of the standard 
argument for n — 1, given e.g. in 0. To prove stationarity of we have to show that 
for all / G C(n L ,R) 



N 



E E E^(^( x ))^(y)(w ;x,x+y )-/^))^^) = - ( 2 - 10 ) 

T]€fl L x,yeA L i=l 

With ()2.6|) it is easy to show that the one point-marginal satisfies 



- /i \ ^(k*, k{ 1) 

= I /un exp^i 



for all i 



<?i(k) 

n, k = (fci, . . . , fc n ) G M n with fcj > 0. For every x, y G Aj, one has 



(2.11) 
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where we introduced the shorthand (771, ?7 n )(x) = (^i(x), ^ n (x)). With this and a 
change of variables in the summation over x we obtain 



N 



ven L i=i x.yeAi, 

9i((Vl,-;Vi+ 1 ,-;Vn)(^-y)) ^ ( {Vl, Vi+ 1, ■■, Vn) (x-y) 



^(^( x -y)) 



Pt(»/(x)) ^ M (r?(x)) 



n P M (f|(*))=0 

zeA L \{x-y,x} 



[2.12) 



which vanishes by (|2.11|) for every / G C(Ql,M.). Thus f|2.10|) is shown and vz is station- 
ary. 

Assume now that v L is an arbitrary stationary product measure of the zero range 
process with generator C. Then 



(£f)„L = for all feC(n L ,\ 



(2.13) 



and, by inserting special functions /, one deduces conditions on the jump rates. Consider 
a configuration f] where there are k = (&!,..., k n ) G N n particles at a fixed site x and 
the rest of the lattice is empty, i.e. fj(y) = <5 y , x k for all y G A L . From the stationarity 
condition (J2.13j) with f = Xfj we obtain for the one-point marginal of v L 



N 



9i (k) 



g n (k) v(k) = J2"(W,k i -l) 



(2.14) 



i=l 



where we set z/(k) = if hi < for some i G {1, . . . , n}. To get (|2.14j) we used 



>(0, ..,0,k h 0, ..,0) = u(O)H - 



1 



AA^(0,..,0,/,0,..,0) 



(2.15) 



as known from the stationary measure of the one component system (J2.4|) . Now let / be 
the indicator function of r7* ;x,x+y , with r) = (<5 yiX k) yeA as above and ki > 0. Using (J2.13j) 
to ()2.15|) we obtain 



Pl {y)v{Q)\g i {k)v{k)-v{k\k l -l) 



Y^ Pj (-y)u{k^,ki- 1,^-1) 



^■(ej + e^) z/(e* + e 3 -) - z/(e,- 



(2.16) 



with the shorthand k = (k* J , fcj, fcj) and e, G M n the unit vector in direction i. ()2.16|) 
holds for all y G and is obviously fulfilled if the two square brackets vanish individually. 
Under the assumption that they do not vanish, one can easily construct a contradiction 
to (|2.14|) . Thus we obtain 



u(k) = u(k i ,k l -l)/g i (k) 



(2.17) 
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for alH G {1, . . . , n}. Applying (j2.17|) twice in different order for arbitrary i ^ j we get 

_ j/Qg, fc, - 1, fcj - 1) _ Kj^, h - 1, fcj - 1) 
lJ ~ ft (k) ^(k*, *, - 1) ^(k^k^-l) ' [ - ' 

and f)2.5j) easily follows. □ 

Remark: If the jump rates are site dependent and satisfy ()2.6|) with potential G x for 
every site x, the measure defined analogous to ()2.8|) is still stationary, since the terms in 
(I2.12|) cancel for each site individually. However, we did not see how to generalize the 
reverse argument to space-dependent rates. 

2.3 Properties of the stationary measures 

The particle density Ri of component i G {1, . . . , n} is translation invariant and given by 
Ri{v) := <^(x)>^ = <^(0)>^ = <9 W logZ(M) > (2.19) 

o 

as a function of the chemical potentials p. For p G Dp the measure has some finite 
exponential moments as (e '^ ))^ = Z(p + 0)/Z(fi) < oo for sufficiently small 6 G R n . 

Therefore R = (R 1} . . . , R n ) : -> D p is well defined with £> p = R(Z)J C (0, oo) n . Due 

to (|2.19|) . logZ(/it) is monotonic increasing on _D M and the compressibility is given by the 
matrix of second derivatives as 

DR(p):= (dp 3 R t {n))=D 2 \ogZ{ii) = [d^. ]og Z(j*j)= (^(0)^(0)^)^,(2.20) 

where (^(0)^(0))^ := (ffc(0)»fc(0)) - (»fc(0)) (fy-(0)) . Thus L>R(p) is symmetric 
and positive definite, because 

a T - (L> 2 logZ(p))a = 0,^(0)) )° >0 (2.21) 

5S=1 

for all aGl n with |a| = 1. Hence the eigenvalues of D 2 \ogZ(fi) are real and positive, 

o 

which ensures that log Z(fi) is strictly convex and R is invertible on D^. Since DR(/i) 

o 

is also continuous, D p is diffeomorphic to D^- We denote the inverse of R by M = 

o 

(Mi, . . . , M n ) : D p — > D M and define the measure ^ := ^M(p)> wmcn is indexed by the 
particle densities p. There exists an a > such that p G _D p for all |p| < a, so is 
defined for small densities. In many cases it is D p = (0, oo) n , e.g. under the assumption 
(|3.3|) or (|3.4j) in Section 3.1. However, there are also cases where there is no stationary 
product measure for large densities. These systems show an interesting condensation 
phenomenon for large p and have been studied for one-component systems in ^3 QE] • I n 
this case the behavior of Z(p) and R(p) at the boundary dDp is of importance. We will 
not discuss this point any further here. 

The stationary current of component i is given by 

Ji(p) = m(pi) (#>„„ = m(pi) exp [M<(p)] (2.22) 



as a function of the particle densities. Here m(pj) = X] y gA L yP*(y) e denotes the 
first moment of the jump probabilities pi, which in general is non-zero and determines 
the direction of the current. The strength is proportional to exp [Mj(p)] and thus a 
monotonic increasing function of the chemical potential. 

The thermodynamic entropy S(p) of the stationary measure is the convex conjugate 
of log Z (//) given by the Legendre transform 

S(p)= sup ( f^p^i-logZip,)) . (2.23) 

With (|2~22"J) we have for all i <E {1, . . . , n} 

d Pi S{p) = M^p) =\og( gi ) Vp . (2.24) 
Therefore we have the following relation for the determinants, denoted by |..|, 

\D 2 S(p)\ = \DM(p)\ = |£>R(M(p))| _1 = |£) 2 logZ(M(p))| _1 > . (2.25) 

Thus S is strictly convex on D p . Note that due to the structure of the stationary measure 
the densities are given as derivatives of the partition function with respect to the chemical 
potentials. This leads to 

d Pj \og(gi) Up = d Pi \og( gj ) Up , (2.26) 

for all i,j 6 {1, . . . , n}, which can be considered as the macroscopic analogue of condition 
(12. 5|) on the jump rates. 



3 Hydrodynamics 

3.1 The hydrodynamic limit 

We will show that under Eulerian scaling t —> t/L, x ^ u = x/L in the limit L — > oo 
the time evolution of the local particle densities p(t, u) is given by the following system 
of conservation laws: 

d 

«9 t p,(t,u) + ^d Ufc jf(p(t,u)) =0, i = l,...,n, (3.1) 

k=l 

where JMp) is the k-th spatial component of the i-th current J,(p) defined in (|2.22j) 

and u = (u\, . . . ,Ud) G A = (M/Z) d is the continuous space variable. To prove the 
convergence, the dynamics of the zero range process has to be well defined in the limit 
L — > oo which is guaranteed by (see j2 E] ) 

sup sup Igiik 1 , kj + 1) — <?i(k)| < oo . (3.2) 

i,je{l,...,n} fceN 

We also need to impose an extra condition on the stationary measure, which is needed 
for the one block estimate (see [H], Chapter 5), as one important part of the convergence 
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proof. There are two alternatives, the first one is to require that the partition function 
Z(p) is finite for all \i G M. n , which is equivalent to the existence of finite exponential 
moments, i.e. 

< exp[0 • r/(0)]>_ = Z(p + 0)/Z{p) < oo for all 6 eW 1 & D^—W 1 . (3.3) 

Note that this implies D p = (0, oo)", avoiding possible problems in case p(t, u) reaches 
the boundary of D p . Alternatively, we can impose sublinearity of the jump rates, i.e. for 
all i G {1, . . . , n} and all b G (0, oo) n there exist aj(b) G 1R such that 

gfj(k) < aj(b) + b • k and lim Z(p) = oo for all p* G dD p . (3.4) 

The second statement is needed to ensure D p = (0, oo)", since sublinearity does not 
rule out S M". Given a solution p(t, u) of (|3.1|) we denote the corresponding local 
equilibrium measure by uh t % , which will be compared to the time dependent distribution 
Trf' of the zero range process, 

v p{t,.) := II u p(^/ l ) and n t = n o S tL (3.5) 

where S t L is the semi-group S t associated to the generator C speeded up by L. The 
proof of the following theorem is an application of Yau's relative entropy method [H], 
which requires some regularity of the solution p(t, u). In general solutions of conservation 
laws develop shocks after a finite time even for smooth initial data (see e.g. [7]). Thus the 
convergence proof is valid only up to the time T of the appearance of the first discontinuity. 



Theorem 2: (Hydrodynamic limit) 

Let p G C 2 ([0,T] x A, [0, oo)") be a solution of JOp for some T G (0, oo) with smooth 
and bounded initial profile p(0, .), satisfying pi(0,u) > p* > 0, i — 1, . . . , n. Under the 
assumptions \3. 2i) and \3. resp. \3.Jj\j let 7Tq be a sequence of probability measures on 
VL L , whose entropy H (tiq | jA q ^) relative to zA Q , is of order o(L d ) . Then 

H{4K { t,.)) = o(L d ) for all t G [0,T] . (3.6) 



Applying the entropy inequality in the standard way (see Chapter 6), Theorem 2 
implies the following 



Corollary: Under the assumption of Theorem 2, for any smooth test function f : A — ► W, 
t G [0, T] and i = 1, . . . , n, the following limit 

r lim Tl E fWL)Vi(t,x) = [ f(u) Pi (t,u)d d u (3.7) 

holds in probability, where rj(t,x) denotes the time t configuration of the zero range pro- 
cess with distribution . 

The proof of Theorem 2 is close to the ones given in [3 |S] and its most important 
steps will be sketched in Section 3.3. 
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3.2 Properties of the limit equation 

Introducing the matrices D3 k (p) := (^J*)^. one can rewrite (j3.1j) in the quasilinear 
form 

d 

d t p(t, u) + J2 D * k {P(t, u )) 9 Uk p{t, u) = . (3.8) 

k=l 

The current is given by 

DJ k (p) = A k (p)DM(p) with A k ] (p)=5 tJ m k (p t )exp[M t (p)} , (3.9) 

where A k is a diagonal matrix and rrik{pi) is the k-ih space component of m(pi). Since 
DlSA(p) = D 2 S(p) is symmetric and positive definite, we can write for all k — 1, . . . , d 

DJ k = (D 2 S)- 1/2 ((D 2 S) 1/2 A k (D 2 S) 1/2 ^(D 2 S) 1/2 . (3.10) 

Thus DJ k is similar to the real symmetric matrix (D 2 S) 1 ^ 2 A k (D 2 S) 1 ^ 2 1 which implies 
that Ylk=i UJ kDJ k is diagonalizable for all uj = (ui, . . . ,uj d ) e M. d with \uj\ = 1 and that 
(j3.1|) is hyperbolic ^H|- The question of strict hyperbolicity, i.e. whether all eigenvalues 
of J2k=i u kDJ k are nondegenerate, cannot be answered in general. It depends on the 
dynamics of the zero range process and one has to check for each system separately. 
With flQH) and (j3~^|) it is easy to see that S(p) satisfies 

n 

J2d Pi S(p)DJUp) = 9 pM M (P)) forallj = l,...,n, k = l,...,d, (3.11) 

i=l 

provided we set F(/x) = Y^=i m i.Pj) f^j ex P\p L j]- (jSHIj) are the dn defining relations 
for entropy entropy-flux pairs of hyperbolic systems El]- For general systems with 
dn > d+1 these equations are overdetermined and the existence of an entropy entropy-flux 
pair is not guaranteed. However, as we just have shown, for a zero range process the Euler 
equation has always a strictly convex entropy S, defined in (|2.23|) . with corresponding flux 
F. 

Systems of conservation laws with entropy are studied in detail in ^H]- in general, by 
transformation to the so-called entropy variables the quasilinear equation (J3.8|) simplifies 
to a symmetric system. In our case these variables are given by the chemical potentials 
/j,(t,u) := M(p(£, u)) and the derivative of the current with respect to p is even diagonal, 

d 

DR(/x(i, u)) a ( /i(t,u) + ^ A^(t, u)) ^./ifi, u) = , (3.12) 

k=l 

where A k (p) := A k (~R(p)) is defined in fUj) . 

To summarize we have shown that the limit equation ()3.1|) is a hyperbolic system with 
globally convex entropy. General results on the existence and uniqueness of solutions of 
such systems are rare, even for the simple form (|3.12j) . So further results for systems with 
open boundaries are based on solid arguments rather than rigorous proofs and presented 
in Section 4. 
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3.3 Proof of Theorem 2 

The proof follows closely the one given in [3], Chapter 6, and jH] and we only sketch the 
main steps. The only part of the proof where the structure of the stationary measure 
for n-component systems enters is ()3.18|h where we use the symmetry of DM(p). Since 
7rf and u^u v are absolutely continuous with respect to each other and with respect to a 
reference invariant measure z/„, a G (0, oo) n , one can define the density 

,uL M ._ K^) _ TT Z(M( a )) " exp^(x)M t (p(t,x/L)] 

dv ~ x A e £ Z ( M (^' x / L )))fi exp^(x)M,H] { - J 

Let Hz(t) := H(Kt\Vp, t n) be the entropy of 7rf relative to v^, t y To establish the estimate 
(13. 6 j) . we will prove a Gronwall type inequality 

< H L (0) + C [ H L (s) ds + o(L d ) (3.14) 
Jo 

with uniform error bound for all t e [0, T]. The entropy production is bounded above by 

d t H L (t) < ! -i- (LC*^(rj) - d t tf(r,)) drf , (3.15) 
Jn L W{il) v ' 

where C* is the adjoint of L in L 2 (z/^). This inequality is proved in Chapter 6, under 
very general conditions covering our case. Using the regularity of M(p(£, .)), the right 
hand side of (|3.15|) can be rewritten as 

n d 
xeA L i=l fe=l 

n d 

- E EE^Ww^i^WftW*)) - jf(p(t,x/L))) 

xeAi i=l fe=l 

n 

(V^fa))" 1 ^*?) = E a * M ^(^' x / L ))(^( x ) -fH&x/L)) ■ ( 3 - 16 ) 



xeA L i=l 



The right hand side of the first line is a telescoping sum and vanishes up to an error 
0(L d ~ l ). Because of the regularity of M(p) a summation by parts permits to replace the 
local variables by their block averages. They are defined as 

^ = 7WTiy E w(y), ^W = ^tw E 9i(v(y)) (3.17) 

1 ' \x-y\<t K ' |x-y|</ 

for i = {1, . . . ,n} and £ G Z + . Using the hyperbolic system (|3.1|) and the symmetry of 
DM we obtain 

n d 

dtM^pfrx/Lj) =-^^d„j}(p{t,x/L))d Vk M j (p(t,x/L)) . (3.18) 

j=l k=l 
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Inserting (EHTf) and (|3~T%J) in (jUSD , 

^ ™ xeA L i=l fe=i 

„ n d 

+ / E EE^ M *(^' X / L )) /f(^(x),p(t,x/L)) d7r t L , (3.19) 

where 

/*(a,b) = J*(a) - J*(b) - VJf(b) • (a-b) . (3.20) 

A bound for the second term on the right hand side of ()3.19|) comes from the entropy 
inequality (see [2], Chapter 6). Dropping the argument of we get 

„ n d 

/ E EE^ M ^> X / L )) ^CH^v^+OiLH- 1 ) . (3.21) 
xeA L i=i k=i 

The first term is estimated integrated in time by using the so-called one block estimate 
(see Chapter 5) 

pt r n d 

lim lim L- d / / E EE (m fc (p0^(»7(x))- ^(^(x))) thf ds = . (3.22) 

This is the only part of the proof where either one of the two regularity assumptions on 
the stationary measure ()3.3|) or ()3.4j) is used. Inserting ()3.21|) and ()3.22|) in ()3.19|) we 
obtain ()3.14|) via ()3.15|) and Theorem 2 follows. 

4 Stationary solutions for systems with open bound- 
aries 

4.1 Uniqueness criterion for the physical solution 

In this section we turn to one of the motivations of our study and apply the results 
of the previous sections to determine stationary density profiles of systems with open 
boundary conditions. General results on existence and uniqueness of solutions to systems 
of hyperbolic conservation laws like ()3.1|) and ()3.12j) are not available and we have no 
other choice than to base our results on solid arguments rather than rigorous proofs. We 
explicitly state our assumptions. 

Although (J3.12)) is derived only up to the first discontinuity, we assume in the following 
the validity of the Euler equation in the sense of weak solutions. They are in general not 
unique and we have to find a criterion to single out the physical solution, i.e. the one 
which describes density profiles of the underlying zero range process. There has been a 
lot of work on this problem (see e.g. [7]) and one possibility is to add a viscosity term 
with a small parameter e on the right hand side of (J3.12)) . A natural choice of this 
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term is the diffusive correction which was neglected in the derivation of (|3.1|) . where the 
small parameter is interpreted as e = (D(L~ d ). In our context this term is given by the 
Green-Kubo formula for the corresponding reversible zero range process with symmetric 
jump probabilities (see Chapter II. 2. 2). Since the symmetric zero range process is a 
gradient system, the viscosity is already determined by the stationary measure itself. We 
get the following dissipative equation 

n d d 

Y DRijilS) dtl*j + Y m k(Pi) 9 u k ex P[/4] = e Y ° kl ^ d ^ui ex P[/4] , (4- 1 ) 

3=1 k=l k,l=l 

for % = 1, . . . , n, with (JkiiPi) = SxeA X fe x /Pi( x )- The solution fi e (t, u) of Equation ([4.1]) 
is unique and, if it converges in a proper sense for e \ 0, the limit /x°(t,u) is a weak 
solution of ([3.12]) . This convergence is in general difficult to show and in the following we 
will just assume that it holds, for details see ^_| Chapter 3.8 and references therein. 

Since D 2 S(p)(5ijaki(pi)d Pj exp [Mj(p)l ) is similar to a symmetric, positive definite 
v _____ ' v 

matrix (analogous to ([3.10]) ). we have for some 5 > and arbitrary a*, e M n , k — 1, . . . , d, 

d d 

Y alD 2 S(p)(6 tJ a M (p t )d P] ex P [Mi(p)]) a* > SYlMl > • (4-2) 

k,i=i %] k=i 

This expresses the viscous dissipation for the entropy and ensures (see [T_j, Chapter 3.8) 
that p°(t,u) satisfies the entropy inequality 

d 

d t S(R(fi°(t,u))) +Ydu k F k {»°( t > u )) < » ( 4 -3) 

k=l 

as to be expected for a physical solution. Thus we expect that the zero viscosity limit 
of ([4.1]) describes the macroscopic chemical potential profiles of the zero range process. 
Note that for stationary solutions the system ([4.1]) decouples and stationary profiles can 
be obtained very easily. They are only determined by the first and second moment of 
the jump probabilities and the boundary conditions, whereas they are independent of the 
jump rates. The rates only enter the partition function Z(fjL) and thus the transformation 
to density profiles via R(/x) ([2.19]) . which is illustrated in the next section. 



4.2 Stationary profiles for one-dimensional systems 

There has been considerable activity to understand the structure of the nonequilibrium 
steady state of systems with open boundaries. Here we study this issue in one dimension 
on the level of the hydrodynamic equations. It turns out that the fugacity variables 
4>i := exp[//j] are the most convenient choice. As a consequence of ([4.1]) . for d = 1 the 
stationary fugacity profiles </>i(u), u G [0,1] are the limit solutions for e \ of the 
equation 

m(pi) d u (t)\{u) = e a{pi) d 2 u (pl{u) , i = 1, . . . , n . (4.4) 
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Figure 1: Stationary solution of (|4.4j) with e = 0.02 for a one-dimensional two-component 
system with open boundary conditions ^i(0) = 1.5, 0i(l) = 1.2, <p 2 (0) = 2, 02(1) = 1-2 and 

jump rates (|4.6|) with c\ = 1.3, C2 = 1. Fugacities: dashed lines (f>\ ( ), (fi 2 ( _ _ ) given by 

(|4.5|) . densities: dash-dotted lines pi ( ), p 2 ( )• Left: m(pi) = 1, m(p2) = 0.5. Right: 

m (Pi) = 1) rn(p2) = —0.5. 

The equations are decoupled and for a system with open boundary conditions <fii(u) = 
exp [Mj(p(w))] for u — 0, 1, the solution is given by 

= 0,(0) + (0,(1) - 0,(0)) ( gi (e)« - l)/(<ft(e) - 1) , (4.5) 

where gj(e) = exp [m(pj)/(ecr(pj))l . The profiles ()4.5|) have a very simple structure. They 
are bounded above and below by </>j(0) resp. </>j(l) and for e \ they converge pointwise 
to flat curves $(u) with a jump at one of the boundaries, if </>j(0) ^ ^ (1) . The location of 
the jump depends on the sign of m(pi), which corresponds to the direction of the current. 
The coupled transformation to the stationary density profile p*(u) involves the fugacities 
of all components and is given by p° = i?j(log05; • • • > log^n) defined in ()2.19|) . 

We illustrate our result with a simple example with two components, for which the 
jump rates are given by the potential 

G(k u k 2 ) = k 1 \ogc 1 + k 2 \ogc 2 + \og(k 1 + k 2 )\ (4.6) 

and thus read gi(ki,k 2 ) = Xfc 4 >oCj(^i + k 2 ), i = 1,2. The two types of particles just 
move with different speeds, but in this case one can calculate the partition function (J2.9J) 
analytically. In terms of the fugacities 0i and cj) 2 it is given by 

Z(<j>i, 4> 2 ) = { <h/ci-<h/c2 ' ri/ r . (4.7) 

[ exp[0i/ci](l + 0i/ci) , for 0i/ci = 2 /c 2 

The transformation to densities is then obtained via Ri((pi,4> 2 ) = 0jC^Z(0i, 2 ) (cf. 
f!2.19jl ). The jump probabilities are chosen such that a{pi) = cr(p 2 ) = 1 and m(pi) = 1. 
Therefore the bulk value <j>1(u), u G (0,1) is equal to the left boundary 0i(O). For 
m(p 2 ) = 0.5 the same is true for 2 and also the density profiles are determined by their 
left boundary value. This can be seen in Figure 1 (left), where we plot the profiles for 
small e = 0.02, for better illustration. For m(p 2 ) = —0.5 the particle species are driven in 
opposite directions, leading to a combination of fugacities in the bulk, which at neither of 
the two boundaries is present. So the bulk densities no longer agree with their boundary 
values, cf. Figure 1 (right). 
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4.3 Concluding remarks to steady-state selection 

The eigenvalues A,(p) of the current derivative DJ(p) defined in ()3.9j) determine the 
characteristic velocities at which small perturbations of a flat fugacity profile propagate. 
For zero range processes the sign of Aj(p) is fixed by the sign of the first moment of 
the jump probabilities m(pi) and independent of the densities p. So no matter how the 
boundary conditions are chosen, the qualitative behavior of the stationary profiles does 
not change. This is in contrast to one-dimensional driven lattice gases with exclusion 
dynamics, where a change of sign of the characteristic velocities is the key ingredient 
for boundary induced phase transitions (see CD])- So our analysis shows that these 
phenomena are not present in zero range processes and the selection of stationary states 
in terms of fugacities is particularly simple. The construction of stationary profiles as 
shown above can be readily generalized to higher space dimensions. Despite the absence 
of boundary phase transitions, the zero range process is a very important interacting 
particle system, last but not least due to its close relation to exclusion models |17j . It is 
one of the few examples of multi-species systems, where the selection of stationary states 
is well understood under very general conditions. 
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